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ABSTRACT 

We investigate the dynamical evolution of compact young clusters (CYCs) near the Galactic center 
(GC) using Fokker-Planck models. CYCs are very young (< 5 Myr), compact (< 1 pc), and only a 
few tens of pc away from the GG, while they appear to be as massive as the smallest Galactic globular 
clusters lO'^Mo). A survey of cluster lifetimes for various initial mass functions, cluster masses, and 
galactocentric radii is presented. Short relaxation times due to the compactness of CYCs, and the strong 
tidal fields near the GG make clusters evaporate fairly quickly. Depending on cluster parameters, mass 
segregation may occur on a time scale shorter than the lifetimes of most massive stars, which accelerates 
the cluster's dynamical evolution even more. When the difference between the upper and lower mass 
boundaries of the initial mass function is large enough, strongly selective ejection of lighter stars makes 
massive stars dominate even in the outer regions of the cluster, so the dynamical evolution of those 
clusters is weakly dependent on the lower mass boundary. The mass bins for Fokker-Planck simulations 
were carefully chosen to properly account for a relatively small number of the most massive stars. We 
find that clusters with a mass < 2 x 10'' M0 evaporate in < 10 Myr. Two CYCs observed near the GG 
— the "Arches cluster" (GO. 121+0. 17) and the "Quintuplet cluster" (AFGL2004) — are interpreted 
in terms of the models; their central densities and apparent ages are consistent with the hypothesis 
that they represent successive stages of cluster evolution along a common track, with both undergoing 
rapid evaporation. A simple calculation based on the total masses in observed CYCs and the lifetimes 
obtained here indicates that the massive CYCs comprise only a fraction of the star formation rate (SFR) 
in the inner bulge estimated from Lyman continuum photons and far-IR observations. This is consistent 
with the observation that many stars in the inner bulge form outside the large clusters. 

Subject headings: celestial mechanics, stellar dynamics — Galaxy: center — methods: 
numerical — galaxies: star clusters 
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1. INTRODUCTION 

The inner few hundred pc of the Galactic bulge 
(the inner bulge) contains stars of a variety of ages, 
in addition to an apparently coeval population of 
galactic-age stars. The evidence for recent massive 
star formation there has grown with the observations 
of several clusters of emission-line stars in that region. 
Intermediate-age stars have also been observed in the 
inner bulge in the form of OH/IR stars (Lindqvist, 
Habing, & Winnberg 1992; Sjouwerman et al. 1998). 
Noting the presence of these young and intermediate- 
age stars, Serabyn & Morris (1995) argued that star 
formation in the central molecular zone occupying the 
inner bulge has been sustained over the lifetime of the 
Galaxy in spite of the inhospitality of that environ- 
ment for star formation, owing to large gas tempera- 
tures and turbulent velocities, strong magnetic fields, 
and strong tidal forces. 

About half a dozen distinct sites of recent star for- 
mation have been found in the inner bulge, most no- 
tably including the "Arches cluster" (G0.12H-0.017; 
Nagata et al. 1995; Cotera et al. 1996; Serabyn, 
Shupe, & Figer 1998) and the "Quintuplet cluster" 
(AFGL2004; Okuda et al. 1990; Nagata et al. 1990; 
Glass, Moneti, & Moorwood 1990; Figer, McLean, & 
Morris 1995, 1999), both lying within ^ 35 pc, in pro- 
jection, of the Galactic center (GC). These clusters 
are very young (< 5 Myr), and compact (< 1 pc), 
while they appear to be as massive as the smallest 
Galactic globular clusters (~ 10* Mq). Their young 
ages are manifested by the presence of ample massive 
stars: the Arches and Quintuplet contain about 120 
and 30 stars having initial mass larger than 20 M0, 
respectively (Serabyn et al. 1998; Figer et al. 1999). 
Morris (1993) suggested that the non-standard star 
formation environment near the GC may lead to an 
initial mass function (IMF) skewed toward relatively 
massive stars (flatter IMF) and having an elevated 
mass cutoff. However, the observational limit for the 
lowest mass in these clusters is still as high as several 
Mq, so not much is known about the lower end of 
their mass function. 

The Arches and Quintuplet clusters are two mem- 
bers of what appears to be a special category of clus- 
ters; indeed, their large masses place them at the 
lower end of the category of super star clusters (e.g.. 
Ho & Filippenko 1996). Another known, potential 
member of this category near the GC is the young 
cluster in the central parsec, but because of its un- 



usual location at the bottom of the Galactic gravita- 
tional potential well, and in the immediate vicinity 
of the central supermassive black hole, it may have 
had quite a different origin. Thus, it is not clear 
that this cluster should be categorized alongside the 
Arches and Quintuplet, even though it has a com- 
parable mass (Figer et al. 1999). In any case, the 
calculations presented in this paper do not apply to 
the central cluster. The stellar cluster in Sgr B2 is 
potentially another supercluster currently in forma- 
tion (e.g., Gaume et al. 1995), although it is not yet 
clear whether it will attain the status of the two more 
evolved clusters. 

With the observations of these compact, young 
clusters (CYCs) near the GC, a natural question 
arises: Why don't we observe older examples of 
CYCs? Are we witnessing very unusual events that 
happened to have recently taken place in the GC, or 
are older clusters rare because they evaporate in a 
relatively short period of time? The proximity of the 
CYCs to the GC implies a strong tidal field, and the 
short relaxation time due to their compactness makes 
very massive stars play an important dynamical role 
through mass segregation, which accelerates the dy- 
namical evolution even more, before these stars disap- 
pear with mass loss and explosions. These two char- 
acteristics suggest fairly short lifetimes for the CYCs. 
The ages of the Arches and Quintuplet clusters are 
estimated to differ by about a factor of two, provid- 
ing an excellent opportunity to compare observations 
with the numerical study of their dynamical evolu- 
tion. Because their rather extreme parameters have 
only recently been recognized, these clusters have not 
yet been numerically and theoretically explored. In 
the present study, we investigate the clusters' life- 
time, iet,, for various parameters plausible for CYCs, 
and compare the results with current observations of 
CYCs and with estimates of the star formation rate 
in the inner bulge. 

A discussion of the time scales related to the dy- 
namical evolution of spherical stellar systems and on 
our choice of simulation methods is given in § ^. We 
describe the details of our models in § |[ and show 
the simulation results in §14. The implications of our 
results are discussed in § |^ and § ^ summarizes our 
findings. 

2. BACKGROUNDS 
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2.1. Time Scales 

The dynamical evolution of stellar systems is af- 
fected by many factors such as galactic tidal fields, 
the stellar mass function, and stellar evolution, as well 
as their overall mass and physical extent. In general, 
a cluster loses its mass by two-body relaxation and 
stellar evolution (via supernovae and stellar winds), 
and the resulting dynamical evolution of the cluster 
is accelerated by the presence of a strong tidal field. 

The evolution of star clusters in a static tidal field 
has been studied extensively. Henon (1961) obtained 
a self-similar solution for an expanding cluster assum- 
ing that there is a central energy source (which he 
speculated to be binaries). Such a cluster is found 
to evaporate completely in 22.5 times the half-mass 
relaxation time, which is defined as 
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where M is the total mass of the cluster, Rh the half- 
mass radius, rfi the mean mass of each star, and G 
the gravitational constant. The In A term appears 
because of the finite size of the cluster, and is called 
the Coulomb logarithm. 

The numerical integration of the Fokker-Planck 
equation by Lee & Ostriker (1987) confirmed Henon's 
general result. These works have been extended to 
multi-mass models for typical globular cluster param- 
eters by Lee, Fahlman, & Richer (1991) and Lee & 
Goodman (1995), and the lifetimes of the clusters in 
units of half-mass relaxation time have been found to 
be shorter than those of single mass models by about 
a significant factor. Thus, clusters with a realistic 
mass function would survive less than 10 tru. 

The evaporation of stars from the cluster is due to 
relaxation, but it takes some time for a star to escape 
from the cluster. Since this "lingering time" is of or- 
der of an orbital period of the star, the evaporation 
of a cluster depends on both relaxation time and dy- 
namical time. The average orbital time scale for stars 
in the cluster can be expressed as: 
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where Rt is the tidal radius of the cluster. 
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Here, Rg is the galactocentric radius, and Mg the 
enclosed galactic mass within Rg. Equations (|l|) and 
(^ give a relation t^yn/trh oc iV^^, where N is the 
number of stars. Thus the evaporation depends more 
sensitively on TV when N is small. Note that since Mg 
is a function of Rg (see eq. |^] ) , equations (||) and 
imply that tdyn can be expressed in terms of Rg only. 

Stellar evolution affects a cluster's dynamical evo- 
lution in several ways. First, mass loss by stellar 
evolution (through stellar winds and explosions) itself 
accelerates the mass loss of a cluster and thus short- 
ens the cluster's lifetime. Furthermore, the mass loss 
by stellar evolution, together with tidal evaporation, 
may prevent the core collapse. Since stars bound to 
a cluster have negative energy and massive stars are 
predominantly located in the core due to mass seg- 
regation, the mass loss from the cluster thus has an 
indirect heating effect. 

The effect of very massive stars may be neglected 
for globular clusters, which have initial tru of ~ 
2^09-10 yj.^ rnuch longer than evolutionary time scales 
of high mass stars. The low mass stars eventually un- 
dergo evolution on a much longer time scale, and their 
effects on dynamics should be very small. Therefore, 
in old globular clusters, the stellar evolution should 
have been effective only in the initial phase when the 
relaxation did not play any role. The subsequent evo- 
lution should have been governed mostly by the re- 
laxation process alone. Thus the study of globular 
cluster evolution may assume that clusters start with 
upper mass boundary, mu, much lower than the real 
initial value (for example, for uriu, Chernoff & Wein- 
berg 1990 adopted 15 with stellar evolution, and 
Lee, Fahlman, & Richer 1991 adopted 0.8 without 
stellar evolution). However, the sizes of CYCs near 
the GC are very small for their masses, leading to ini- 
tial trh as short as ^ 10^^^ yr. Moreover, the mass 
segregation time, trhrh/mu, may be even shorter than 
the lifetime of the most massive star. Then very mas- 
sive stars play an important role in the evolution of 
the cluster through rapid mass segregation and sub- 
sequent demise. 

To summarize, the dynamical evolution of a clus- 
ter is determined by relative differences between time 
scales for the above processes, two-body relaxation in- 
cluding core collapse and subsequent expansion, tidal 
evaporation, stellar evolution, and mass segregation. 
For example, stellar evolution determines the cluster 
evolution time scale when the chemical evolution of 
massive stars occurs much faster than the relaxation 
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time scale of the cluster, and the relaxation governs 
the cluster evolution when its time scale is shorter 
than that of stellar evolution. Also, if a cluster has 
a large stellar mass range, the segregation of massive 
stars toward the central parts shortens the relaxation 
time in the core and makes relaxation processes more 
important. Finally, a strong tidal field accelerates 
all these processes. Since the time scales of these 
processes scale with cluster parameters in different 
ways, scaling the evolutionary route of one cluster to 
another is quite limited, and numerical simulations 
are generally required for the unexplored parameter 
regime. 

2.2. N-body vs. Fokker-Planck Models 

For the study of dynamical evolution of spheri- 
cal stellar systems, N-body and Fokker-Planck sim- 
ulations are among the most widely used methods. 
The N-body method is certainly more realistic and 
requires less assumptions than other methods, but 
real-number N-body simulations for CYCs having 
N ~ 10^~^ are still prohibitively expensive in terms of 
computing time, especially in the context of a simu- 
lation survey in which a range of parameters is be- 
ing explored. Many studies based on the N-body 
method use a time-scaling technique to simulate a 
large-number system with a smaller number of stars 
on behalf of computing cost. However, correct scal- 
ing is not always guaranteed when the stellar evolu- 
tion time scale, tse, is not well separated from tdyn 
and trh, as might be the case of CYCs (see, however, 
Aarseth & Heggie 1998 for a compromise technique, 
"variable scaling"). Also, as mentioned above, the 
evaporation rate of tidally-limited clusters is a func- 
tion of both trh and tdyn , which scale differently with 
N and M. Thus, direct N-body calculations of the 
system with a smaller number of stars are not able to 
correctly mimic the evaporation of the system with a 
larger N. Moreover, when the mass function is steep, 
the upper mass boundary, r7i„, can be a function of 
the total number of stars in the cluster for statisti- 
cal reasons (the more stars, the higher m„). Then 
the system with a smaller number could not properly 
represent the full mass range of the larger system. 

On the other hand, although the Fokker-Planck 
method requires more assumptions and approxima- 
tions, it is much less expensive and gives statistically 
correct results. These features are critically benefi- 
cial, not only when the study is meant for a parameter 
survey, but also when statistically stable results with- 



out random noise are necessary for theoretical anal- 
yses. For these reasons, as a first step of our study 
of the fate of CYCs, we simulate the dynamical evo- 
lution of CYCs using Fokker-Planck models. We are 
planning N-body simulations for a few representative 
cluster parameter sets in a foUowup study. 

A recent comparison between Fokker-Planck and 
N-body simulations for tidally-limited clusters showed 
a discrepancy in cluster lifetime between the two 
methods (Fukushige & Heggie 1995; Portegies Zwart 
et al. 1998). However, Takahashi & Portegies Zwart 
(1998) were able to show that this discrepancy can be 
successfully removed by adopting anisotropic Fokker- 
Planck models with the "apocenter criterion" and an 
appropriate constant for the speed of star removal 
behind the tidal radius. Here we adopt anisotropic 
Fokker-Planck models and the apocenter criterion as 
an effort to minimize the possible discrepancy with 
N-body simulations. 

3. MODELS 

We have used the two-dimensional (energy-angular 
momentum space), orbit-averaged Fokker-Planck code 
by Takahashi (1997). The code assumes that star 
clusters are spherically symmetric and are in dynam- 
ical equilibrium. 

For our initial models, we adopt a multi-mass dis- 
tribution function composed of single mass King mod- 
els with equal velocity dispersions. It is assumed that 
clusters initially fill the tidal radius and that all stars 
in the cluster are formed simultaneously ait = 0. The 
relation between M and Rt is determined via equa- 
tion (||) and 

which was adopted from Genzel & Townes (1987). 
Primordial binaries are not considered in this study, 
and the Coulomb logarithm is taken as In A = In N. 

Heating by binaries formed through three-body 
processes and tidal-capture processes is included. For 
the heating rate per unit volume by three-body bina- 
ries, we adopt the formulation by Lee et al. (1991): 

E,,^A.21xlO^G'(^J2^) -c, (5) 

where the summation is over all mass components, 
rii and vf are the number density and velocity dis- 
persion of component i, respectively, and is the 
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mass-weighted, central velocity dispersion. On the 
other hand, we use a modified version of the tidal- 
capture binary heating rate per unit volume by Lee 
& Ostriker (1993): 

Etc = ^ + 'rnc)o-tcVtjninj(j)c, (6) 

i j>i 

where the summation i is only for main-sequence 
stars, ffic is the mean mass in the center, ate the 
tidal-capture cross section, Vij the relative rms veloc- 
ity between components i and j, and (pc the central 
gravitational potential. The fitting formulae of Kim 
& Lee (1999) were adopted for ate- 

The effect of the Galactic tidal field is assumed to 
be constant (circular orbit around the GC and spher- 
ically symmetric potential due to the tidal force) and 
is realized by imposing a tidal boundary. A star is re- 
moved if its apocenter distance, Ra, exceeds the tidal 
radius, Rt (apocenter criterion; Takahashi, Lee, & In- 
agaki 1997). Removal of stars satisfying the apocenter 
criterion is realized by the formalism of Lec & Ostriker 
(1987): 



dt 



1- ^ 



-I 1/2 



(7) 



where / is the distribution function, Et the tidal en- 
ergy, and pt the mean mass density within the tidal 
radius. This formula accounts for the persistence time 
of an escaping star, and the speed of star removal is 
determined by a dimensionless constant aesc- Fol- 
lowing Takahashi & Portegies Zwart (1998), we have 
adopted aesc = 2. 

The effect of stellar evolution is realized by sec- 
ularly decreasing the mass of each component. Fol- 
lowing Chernoff & Weinberg (1990), the mass of each 
component is linearly decreased by 



mt{t) = mi(0) - [TOi(O) - ruij] 



for 



< t < T. _ 



(8) 



where mi{t) is the mass of component i at time t, rriij 
the final mass of the component, t^^i the lifetime of 
a star with mass ^i±^ (see below for ?Tjj^i). For 
t < T^_^_i , mi{t) has a constant value of mi(0), and for 
t > Tj_i, mi{i) is set to rn;./. For intermediate and 
low mass stars, the mass evolution between the main- 
sequence and degenerate stage is so abrupt that it can 



be approximated with a step function. However, in 
case of massive stars (m > 20 Mq), the mass evolu- 
tion is continuous throughout the stellar lifetime, so 
it is quite difficult to define the time for which a star 
maintains its initial mass and after which the star 
becomes degenerate. Furthermore, it is not trivial 
to implement such continuous mass loss into a mass 
component which is a representation of a set of stars 
with a certain range of masses and mass loss rates. 
For this reason, here we define an "effective lifetime" 
of a star such that 

rt, 



1 



•m[t) dt, 



(9) 
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where •m[t) is the mass evolution of a star with initial 
mass m(0), tf is the time required for the mass of a 
star with initial mass m(0) to become the final mass 
mf. We have adopted Schaller et al. (1992) for m{t) 
and i/, and Drukier (1995) for m/. Figure]^ is a plot 
of Teff as a function of m(0) obtained in this way. 
Then this Tf,ff is used for Tj^i in equation (^). 

The mass spectrum is represented by a set of dis- 
crete mass components. For the mass binning, equal 
logarithmic space binning has been used in many pre- 
vious studies. However, the same logarithmic binning 
is not appropriate for representing largely different 
mass spectrum slopes at the same time. In particu- 
lar, the logarithmic binning may result in too small 
number of stars for the most massive bins of a clus- 
ter with steep mass spectrum and/or with small Af , 
which is not desirable for Fokker-Planck models. One 
could lower the upper mass boundary in accordance 
with the cluster parameters to avoid such a problem, 
but then the problem becomes how one "decides" the 
upper mass boundary. Generally, one is interested in 
a cluster with given M, a, and m/. In cases where 
the relaxation time of a cluster is much longer than 
the stellar evolution time for massive stars, the effect 
of those massive stars on the dynamical evolution of 
the cluster may be neglected and the choice of ruu is 
relatively unimportant. However, in case of our tar- 
get clusters, the role of massive stars is expected to 
be important because of their short relaxation times, 
and thus rUu should be chosen carefully. For this rea- 
son, we bin the mass in a way that guarantees all mass 
components to have the number of stars larger than 
a certain minimum value. The mass range is binned 
to make all components have equal value of 



N{m) dm, 



(10) 
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where [m^_i^m^j^i \ is the mass range of a component 
i. Thus the mass ranges of each component differ for 
models with different parameters. For a given M, rUu, 
mi, a, and iVf,i„ (number of bins), one can determine 
the value of (3 that assigns a stellar number of Nu to 
the most massive bin. Once (3 is determined, m^^i 
is obtained by equation ( |l0[ ) with given to/ and m„. 
Here we adopted m„ = 150 Mq, mi = 0.1, IM©, 
.^bira — 15, and = 50. In this way, the effects 
of the most massive stars can be realized in Fokker- 
Planck models correctly and consistently. Note that 
/3 = 1 is for equal mass binning, and /3 = is for equal 
number binning. 

We choose the initial mass spectrum to be a simple 
power law 

dN[m) (xm-"^ dm, (11) 

where dN{m) is the initial number of stars with 
masses between to and to + dm, and a = 2.35 gives 
a Salpeter initial mass function. The initial mass of 
component i, TOi(O), is calculated by 



m,(0) 



J m . 1 
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dm 



(12) 



and the initial total mass of component i is 



M,(0) 



m 1 
'+3 



m . 1 
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N(m)m dm 



(13) 



The number of stars in component i, Nt, is then de- 
fined to be Mi/ mi. In Table 0, /3 and the initial mass 
of the most massive component, toi5(0), are given for 
all models. Note that (3 and tois are determined by 
M, TO/, and a because we fix to„, Ni,in, and N^- 

4. RESULTS 

An example of the cluster evolution in our mod- 
els is shown in Figure ^. A mild core collapse takes 
place at ~ 0.6 Myr, and the postcoUapse expansion 
is accelerated by indirect heating due to stellar evo- 
lution starting from ^ 2.2 Myr, which is clearly man- 
ifested by a rather abrupt change in Rt slope. The 
mildness of the core collapse is only an apparent phe- 
nomenon: the relative increase of pc during the col- 
lapse becomes larger when pc is plotted for only the 
most massive component. While the core collapse and 
postcollapse expansion time scales are determined by 
the size, structure, and mass function of the cluster. 



the stellar evolution time scale is dependent only on 
the mass function. Since these time scales are compa- 
rable for most of our models, the lifetime of a cluster 
will be determined jointly by these time scales in a 
complex way. In this section, we discuss the effect of 
each cluster parameter on the evolution of the clus- 
ter. We pay special attention to t^v, which is given in 
Table @, since this is closely related to the fate of the 
cluster. 

Most of the discussion of evaporation and dynam- 
ical evolution assumes that the time is expressed in 
units of relaxation time. That is useful if the rele- 
vant process is related mostly to the two-body relax- 
ation time. However, the stellar evolution is impor- 
tant in the situation that we are considering here, as 
discussed in § 2^ . Therefore, it is more convenient to 
express the results in absolute time for the CYCs. 

4.1. Initial Concentration {Wq) 

King models, which we use as initial models, are 
a one-parameter family, whose parameter deter- 
mines the initial degree of central concentration (the 
larger the value of Wq, the higher the concentration). 
Figure H shows the evolution of models 101, 111, and 
112, which have the same initial conditions except 
for Wq. While the times to total evaporation tf,y are 
almost the same, the epochs of the central density 
peak for these three models are largely different. The 
time to core collapse, t^c, is determined by the central 
relaxation time, trc, which is smaller for more con- 
centrated clusters, while the tidal evaporation rate is 
determined by a cluster's global properties. Clusters 
slowly expand in the postcollapse phase, and this ex- 
pansion accelerates the tidal evaporation. The evap- 
oration times are nearly the same for these models 
because we required our models to have the same tdyn 
(the same Rg), M, and Rt. The evaporation is not 
sensitive to the detailed structure of the cluster, but 
depends mainly on macroscopic parameters such as 
M and Rt- Thus the insensitivity of tev on Wq will 
hold for the parameter regime covered by our other 
models as well. The evolution of the central param- 
eters depends on Wq just because the core collapse 
is determined by the local conditions in the center, 
which depend on the initial models. 

4.2. Mass Range (to/) 

The effects of lower and upper mass boundaries on 
tev might be seen by comparing models 101, 102, and 
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103. One may expect that mi is more important in 
determining t^v than m„ because fh varies more with 
mi than for our a values. However, interestingly, 
the results of models 101 to 103 show the opposite 
behavior: t^v is approximately proportional to m^^, 
but is almost insensitive to mi (see next subsection 
for the relation between tev and mis). Generally, a 
cluster with a mass spectrum rearranges itself so that 
heavies dominate the central region of the cluster and 
lights the outer region in a time scale ~ trhmi/mu 
(~ trhmi/miQ in our case; mi5 better represents the 
upper mass boundary than to^ does because mis ac- 
counts for the mass spectrum; see §P). However, the 
mass ranges of the above models are so wide (thus 
the ratios mis/m; are large) that even the outer re- 
gion is dominated by heavies. This can be clearly 
seen in the density profile plot (Figure ||), where heav- 
ier components are dominant throughout most of the 
cluster. This presumably results from the selective 
evaporation of lighter components due to equiparti- 
tion. However, the dominance of heavies in the outer 
region does not solely explain the similar t,,^ of models 
101 and 102 because the masses of the light compo- 
nents of these models are largely different and may 
contribute to trh in dissimilar ways. 

The postcoUapse density profiles are not far from 
the profile at the end of the collapse, and the latter 
may be described by the empirical relations: 



pM cx B(3/2,p, + l)(-£;o)P' 
din Pi 



d\nr 

Pi 



-0.23(p, -t- 3/2) 



Pu 

mu 



^3/2 (14a) 
(14b) 
(14c) 



where subscripts i and u are for the ith compo- 
nent and the most massive component, respectively, 
B(x,y) is the Beta function, and Eq the central po- 
tential (Bahcall & Wolf 1977, Cohn 1985). The above 
relation stems from fi oc {—E)P\ which is an approx- 
imate distribution function that a cluster has dur- 
ing the collapse. Cohn fomid that pu « 8.2. Equa- 



tions (14a) and (14b) give nearly the same profiles 
for Pi ^ 1, which is seen also in Figure |^. Thus when 
a cluster has a wide mass range, its outer region is 
dominated by heavier components, and t^h is not de- 
pendent on mi as long as mi <^ m„/8.2 (<^ mi5/8.2 
in our case). Such insensitivity of trh to mi is also 
seen between models 113 & 115 (smaller a), models 
111 & 117 (smaller Wo), models 112 & 118 (larger 
Wo), models 125 & 129 (smaller M), and models 122 



& 127 (larger M). On the other hand, t^v's of mod- 
els with a > 2.35 do depend on m/, especially for 
M = 5 X IQ-^ M0 (see Figure ^). This is because their 
mass ranges [mi, mis] are not wide enough. 

4.3. Initial Mass Function (a) 

We discuss the effect of a on tev in terms of mis 
(recall that mis is a function of a). Models 113, 101, 
142, and 114 have the same M (= 2 x 10^ Mg) but 
different a values. As noted above, t^v is inversely 
proportional to mis (see Figure ||) because of the pre- 
dominance of heavier components in determining t^h 
after rapid mass segregation. Models 124, 125, 141, 
and 126 (M = 5 x lO'^M©) show the same propor- 
tionality too. 

Models 121, 122, 143, and 123 (M = 10^ M©; dif- 
ferent a) do not exactly follow the t^y oc m^^ rela- 
tion (although there is still a clear inverse relation 
between tev and tois). As a cluster's mass becomes 
larger, its lifetime and tois increase, and therefore the 
stellar evolution timescale is a relatively smaller frac- 
tion of the cluster's lifetime, and the relative amount 
of mass that is subject to loss by stellar evolution 
before cluster evaporation becomes larger. Since the 
relative amount of mass that is subject to loss by stel- 
lar evolution in a given period is a function of a, tev 
is not simply proportional to for the more mas- 
sive clusters in which the stellar evolution effect is not 
negligible. 

4.4. Galactocentric Radius {Rg) 

Without stellar evolution, clusters with the same 
M would evolve with a time scale proportional to 
trh oc R^/^ oc (see eqs. §, g, and g). How- 

ever, the effect of stellar evolution on tev is important 
in many of our models, and the applicability of this 
relation is limited. Figure ^ shows tev of models with 
a = 2 and 2.35 as a function of Rg and M. In case 

of a = 2, the tev « relation holds only for the 

models with M = 5 x 10'^ Mq, where models end be- 
fore any significant stellar evolution, due to relatively 
low mis values. Stellar evolution significantly affects 
tev of the other a — 2 models, for which the expo- 
nent in the above relation is smaller than 9/10. For 
a — 2.35 models, where the role of massive stars is 
relatively less important, the effect of stellar evolution 
is much less apparent. Models with Af < 2 x 10^ M© 
end before stellar evolution has had a significant ef- 
fect, and the slope of log tev over logi?g of models 
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with M = 10^ Mq is affected by stellar evolution only 
slightly. 

5. DISCUSSION 

5.1. Overview of the Results 

Our tev survey results show that for our parame- 
ter regime, CYCs will evaporate in less than or about 
10 Myr, except for some cases with M = 10^ Mq and 
a > 2.35. The compactness of CYCs, strong tidal 
fields, and large stellar mass ranges are jointly re- 
sponsible for the short lifetimes of CYCs, but among 
these, the strong tidal fields play the most important 
role. Figure compares the evolution of model 142 
with a model which is the same in every respect, ex- 
cept that Rt is ten times larger (because the density 
profile is the same as in model 142, the cluster does 
not initially fill up the tidal radius in this case). Here, 
the effect of the tidal field is dramatically illustrated: 
when CYCs are located in tidal fields much weaker 
than in the GC, their lifetimes are several orders of 
magnitude longer. This could be the case for the R136 
cluster in the LMC, which has cluster parameters sim- 
ilar to those of CYCs, but is apparently subjected to 
a much weaker tidal field. 

5.2. Comparison vifith Observations 

Of the Arches and Quintuplet clusters, the former 
is significantly more compact, and is estimated to be 
younger and more massive (see Table H) . The smaller 
central density of the Quintuplet may be due to an 
initial condition with lower concentration or smaller 
M than that of the Arches. However, because their 
ages are probably quite different, it is interesting to 
suppose that the two clusters had the same initial 
conditions and to see if they represent two distinct 
epochs in the evolutionary track of a cluster. Since 
the Arches is estimated to have M ~ 1-2 x 10^ Mq 
and Rg ~ 30 pc, our models with M = 2 x 10^ Mq 
and Rg = 30 pc can be considered to represent these 
clusters. Figure ^ shows the evolution of model 142, 
whose parameters are M = 2 x 10^ Mq, a = 2.35, 
Wo = 4:, Rg = 30, and mi = IMq. By comparison 
with the locations of the clusters in this diagram, we 
infer that the Arches is located near the epoch of core 
collapse, and the Quintuplet is near the end of the 
evolution. Figure || lends plausibility to the idea that 
the two clusters started with the same initial condi- 
tions and that their largely different central densities 
are simply due to the age difference. 



The two clusters may equally well be represented 
as different evolutionary stages of other model clus- 
ters having different values of a and m;. Model 114, 
which has the same parameters as model 142 except 
for a = 2.5, gives an evolutionary route similar to that 
of model 142, with a slightly longer t^v Meanwhile, 
smaller a values result in shorter t^v, and smaller mi 
values give longer tev. Thus clusters with an appro- 
priate combination of smaller a and mi will also be 
able to describe the two clusters together. 

The mass function evolves due to mass segregation, 
selective ejection of lighter stars, and stellar evolution. 
Furthermore, the mass function varies with distance 
from the cluster center. Thus observations of the mass 
functions at several different radii, along with a de- 
termination of the cluster age, would all be needed 
to infer the cluster IMF. To illustrate this point, we 
present in Figure the evolution of a for the whole 
cluster and for four equally-spaced annuli of model 
142 . Considering the observational inaccessibility of 
low mass stars, and the rapid evolution of massive 
stars, the determination of a was restricted to the 
mass range [3, 30] Mq. While the whole cluster mass 
function evolves relatively slowly until before the fi- 
nal disintegration phase, the mass function of each 
annulus evolves rather rapidly from the beginning. 
The differences in a values between annuli are largest 
near the core collapse. The a values of the outer an- 
nuli initially grow due to mass segregation, and later 
decrease as the tidal radius shrinks. Mass functions 
of CYCs measured with HST will soon be reported 
(Figer et al. 1999), and the radial dependence must 
clearly be considered. Since the mass function evolves 
from the very beginning, even the current mass func- 
tion of the younger CYC, the Arches, will somewhat 
differ from the IMF. 

The final disintegration of a cluster takes place in a 
relatively short period of time compared to the clus- 
ter's lifetime. Fukushige & Heggie (1995) argue that 
a cluster finally disrupts by losing equilibrium when 
Rh/ Rt exceeds a certain critical value. At this point, 
the central density and velocity dispersion decrease 
abruptly, as may seen in Figure ^. We find that our 
models have Vc less than or about 5kms~^ only dur- 
ing the final disruption phase, before which clusters 
maintain v^s larger than or about lOkms^^. There- 
fore, an observation of Vc for the Quintuplet will be 
able to provide important information on the cluster's 
position in its evolutionary track. 

If the Quintuplet is indeed in a disruption phase, 
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tey values calculated in this study may give a con- 
straint on the initial condition of the cluster: the 
Quintuplet is not likely to initially have a > 2.5 if 
its mi is as low as 0.1 M©, because the Quintuplet 
would not otherwise be in a disruption phase at its 
supposed current age. 

5.3. The Star Formation Rate 

The expected number of observed clusters at a cer- 
tain epoch in the inner bulge, Nobs, is basically a func- 
tion of M, tey, and the star formation rate (SFR) in 
the region: 



Nobs - 



SFR 
M 



Rg)^{M,Rg)dRgdM, 



(15) 

where (j){M, Rg) is a normalized probability that a 
cluster is formed with M at Rg. We assume initially 
that stars are all born in clusters. As shown in Figure 
^ tev should apparently be dependent on both M and 
Rg. Thus here we adopt a simple, order-of-magnitude 
analysis for the relation between t^v and Nobs- First 
we split the variables in t^v and </>. For the former, 
we adopt a simple form. 



Rn 



30 pc 



0.75 



M 



2 X 104 Mq 



(16) 



The virtue of this approximation is that it, along with 
variable splitting for (j), makes equation (15) inde- 
pendent of M, whose distribution is not well known 
(in fact, the behaviour of t^v's is better described 
with cx M^/^ rather than oc M, but we choose the 
latter for simpler estimation). For the latter, we 
define cf) = (j)M{M)(f)R{Rg), where (j)M and 4)r are 
normalized for their variables, and we assume that 
follows the density profile in the inner bulge pg: 
(pR OC RgPg (x const {pg (X R~^ is adopted instead of 
Pg OC Rg^'^ for simphcity). We take [0, 100] pc for 
the range of Rg (the final result is not so sensitive 
to the choice of this range) . Then equation (|l^) now 
becomes 



.r f SFR 
VO.lMoyr ^ 



f-lOOpc 



/Opc 



V30 pc 



0.75 



iMRg)dRg j (j}M{M)dM 



20 



SFR 



(17) 



^O.IMq yr-\ 

Nobs is then equal to the actual current number of 
large, distinguishable clusters in the inner bulge, two. 



if SFR is as low as 0.01 Moyr"^, a value which hap- 
pens to be very close to the mass of the two clusters 
divided by their lifetime, which is a crude estimate of 
the SFR without any consideration of dependencies 
on Rg and M. 

Giistcn (1989) estimated the SFR in the inner 
bulge to be 0.3-0.6 M© yr~^ from the global produc- 
tion rate of Lyman continumn photons, and 0.05 Mq yr~^ 
from the luminosity of the discrete far-IR sources 
measured by Odenwald & Fazio (1984). Glisten's es- 
timates are based on the assumption of a Salpeter 
mass function with mi of 0.1 Mq. For a flatter mass 
function with elevated mi, his SFR becomes smaller, 
although still considerably larger than 0.01 Mq yr~^. 
It thus appears that a considerable fraction of stars in 
the inner bulge forms outside of clusters, or in much 
smaller clusters than the Arches and Quintuplet. 

A number of other sites of star formation are 
known in the GC, but none apparently has a stel- 
lar luminosity within an order of magnitude of the 
Arches and Quintuplet. Partial surveys have been 
carried out at near-IR wavelengths to identify other 
potential clusters near the GC (Gatchpole, Whitelock, 
& Glass 1990; Figer 1995; Philipp et al. 1999), but 
no others have yet been identified, and it seems likely 
that if there are any others comparable to the Arches 
or Quintuplet, they are highly obscured, even at 2 
microns. 

We therefore speculate that the circumstances of 
the formation of these clusters was peculiar in some 
fundamental way, such as catastrophic formation by 
extremely strong shocks, or direct collision of two 
dense molecular clouds. If so, then young stars near 
the GG are formed in multiple modes, and one cannot 
ascribe the total star formation rate there to massive 
clusters alone. 

5.4. Binary Heating and Close Encounters 

We find that in all of our models, heating by three- 
body binaries exceeds that by tidal-capture binaries 
during the core collapse, and that the postcoUapse 
expansion is driven by three-body binaries. More- 
over, in most models, the cumulative number of tidal 
binaries formed is less than unity until the end of evo- 
lution. These numbers are one to two orders of mag- 
nitude smaller than the number of coUisional mergers 
between stars during the first few Myr calculated by 
Portegies Zwart et al. (1999) for the R136 cluster 
in the LMG with initial central mass density sim- 
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ilar to our models. The discrepancy is due to the 
fact that we did not consider the evolution of massive 
stars into a giant phase in which the stellar radii are 
largely increased. Thus our tidal-capture binary for- 
mation rates are somewhat underestimated, but since 
the heating per binary is much larger for three-body 
binaries, the inclusion of radius evolution into our 
calculations would not change the overall dynamical 
evolution of the clusters except only for some cases 
with considerably large initial central densities and 
steep IMFs (these conditions are more favorable to 
tidal-capture binary formation; see, e.g., Kim, Lee, & 
Goodman 1998). In these eases, the postcollapse ex- 
pansion will be driven by tidal-capture binaries and 
clusters will expand more slowly than the ones driven 
by three-body binaries. 

5.5. Dynamical Friction 

Clusters are subject to dynamical friction as they 
go through field stars while orbiting around the galaxy 
center. The drag by dynamical friction causes the 
clusters to lose energy and spiral in toward the center, 
where the tidal field is stronger. When the field stars 
have a density distribution oc Rg^, the time required 
for a cluster initially on a circular orbit to reach the 
center is given by 

^ 2xl09yr ( Rg V ( M 

InA V30pcy V2xlO4M0y 

where Vc is the circular orbital velocity of the cluster 

and InA is the Coulomb logarithm for the field stars 
(Binney & Tremaine 1987). Reasonable values for the 
inner bulge are Vc — 170kms~^ and InA « 10. Then 
tfric would be comparable to or smaller than tev cal- 
culated in this study only for clusters with extremely 
small Rg and/or large M. Thus, most CYCs will 
disintegrate through tidal evaporation before a sig- 
nificant orbital decay due to dynamical friction takes 
place. 

5.6. Limitations 

We assumed that clusters consist of coevally formed 
stars and start their main-sequence stages at the be- 
ginning of simulations. However, stars may gravita- 
tionally interact with others already from pre-main- 
sequence or even from accreting protostar stages and 
these stages may be important to CYCs because of 



their short tgr's- For very massive stars, the time scale 
of these stages prior to the main-sequence must be 
much shorter than tev's calculated in this study, but 
it becomes larger than 10 Myr for stars lighter than 
2M0 (Bernasconi 1996). Thus if lighter stars system- 
atically form earlier than heavier stars, the dynami- 
cal evolution of lighter stars until the formation of 
heavier stars may become important. The role of the 
stages prior to main-sequence on the cluster's lifetime 
is not clear though, because not much is yet known 
about the details of cluster formation, such as what 
fraction of the gas goes into star formation (star for- 
mation efficiency), how quickly the gas left over from 
star formation disappears from the cluster by stellar 
wind and supernova explosions (timescale of residual 
gas expulsion), and which gas clumps (as a function 
of mass) contract into stars first and act as discrete 
gravity sources. These uncertainties are particularly 
large for the unusual clusters we are considering here. 
These issues are in the arena of N-body simulations 
and the first two issues were considered by Tenorio- 
Tagle ct al. (1986) and Goodwin (1997a, b) for yoimg 
globular clusters. No observational evidence for the 
presence of residual gas in the Arches and Quintuplet 
has been found yet. Proper inclusion of these effects 
into early cluster dynamics will be possible only after 
a significant advance in the study of star formation is 
achieved. 

Throughout the calculations, we have set A''i5 = 
50, but this value is arbitrarily chosen. Higher val- 
ues would give more statistically correct results, but 
would not represent the whole mass range well. This 
is one of the intrinsic limitations with Fokker-Planck 
simulations. To sec the dependence of tev on N15, 
we performed models 142 and 145 with iVi5 = 150, 
and obtained only ^ 25 % larger tev values. This im- 
plies that the choice of A''i5 does not significantly alter 
our general results. On the other hand, our single- 
power-law IMFs may be too simple to represent re- 
alistic IMFs. Thus we tried a Kroupa IMF (Kroupa, 
Tout, & Gilmore 1993), which is a three-part power- 
law with a = 2.7 for m > IMq, 2.2 for 0.5 Mq < 
m < IM0, and 1.3 for 0.08 Mq < m < 0.5 M©, for 
M = 2 X IO^Mq, = 30 pc, and Wq = 4. We found 
that the model with a Kroupa IMF has results sim- 
ilar to those of model 116, in wliieli the parameters 
are a = 2.5 and mi = 0.1 Mq. Since the Kroupa 
IMF was derived for the Galactic disk, it may not be 
applicable to the GC where the star formation envi- 
ronment is largely different from that in the disk, but 
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wc believe that even the evolution of a cluster with 
realistic GC IMF can be well approximated with one 
of our simple power-law IMFs. 

Orbit-averages uHc;d in our Fokker-Planck simu- 
lations are valid when tdyn trh- Since trh ~ 
{O.IN/ In A)tdyn, CYCs, which have of ~ lO^-^, 
apparently satisfy the requirement. However, the 
wide mass range of CYCs shortens the relaxation 
time in the core through mass segregation, and the 
above requirement may not hold in the core for some 
of our models, especially the ones with the largest 
mass ranges. In addition, while the Fokker-Planck 
method also assumes tdyn ^ tse, some of our models 
with large M and small a may have periods when the 
assumption is temporarily violated. The validity of 
using the Fokker-Planck method for these marginal 
situations can be assessed using N-body simulations. 

6. SUMMARY 

We have investigated the dynamical evolution of 

CYCs near the Galactic center with anisotropic Fokker- 
Planck models. Stellar evolution and heating by 
three-body and tidal-capture binaries were included 
in the calculation, and the apocenter criterion with 
an appropriate removal speed constant was adopted 
for the tidal evaporation. Mass bins were chosen de- 
pending on cluster parameters to properly account for 
a relatively small number of most massive stars. 

The evolutionary time scales of CYCs are fairly 
short because of their compactness and strong tidal 
fields. For our parameter regime (1.5<a<2.5,m; = 
0.1 & IMq, 5 X IO^Mq < M < IO^Mq, and lOpc < 
Rg < 100 pc), clusters evaporate in < 10 Myr except 
for few clusters with M = 10^ M©. Core collapse 
takes place in most clusters, and three-body binaries 
dominate the heating. 

Unlike globular clusters, very massive stars in 
CYCs play important roles in dynamics because re- 
laxation and mass segregation times are comparable 
to or even smaller than those stars' lifetimes. Rapid 
mass segregation due to a large difference in upper 
and lower mass boundaries accelerates a cluster's evo- 
lution even more. Strongly selective ejection of lighter 
stars makes massive stars dominate even in the outer 
regions of the cluster and the evolution of such clus- 
ters is weakly dependent on the lower mass boundary. 

We found that the Arches and Quintuplet clusters, 
which have quite different central densities, central 
concentrations, and ages, may be placed on one evo- 



lutionary track. Among our cluster models, that with 
M = 2x 10" Mq, Rg = 30 pc, a = 2.35 or 2.5, and 
mi = 1 Mq gives the evolutionary track that best de- 
scribes both clusters simultaneously. The first two 
parameters are from observational constraints, and 
an appropriate combination of smaller a and mi can 
also describe the two clusters together. 

The total mass in observed CYCs is several 10"^ Mq 
and the lifetime of CYCs with observed M at ob- 
served Rg is a few Myr. These numbers result in an 
inner bulge SFR of ~ 0.01 Mq yr~^. A more detailed 
(but still approximate) estimate of the SFR with a 
consideration of distributions in M and Rg happens 
to be similar to this simple estimate. Lyman con- 
tinuum photons and far-IR observations in the inner 
bulge region give SFRs of 0.05-0.5 M© yr~-^. We sug- 
gest that the mode of formation of CYCs is pecu- 
liar, likely resulting from strong shocks or collisions 
of dense molecular clouds. 
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Table 1 

Simulation Parameters and Eesults 



Models 


a 


Wo 


M 
(Mo) 


Rn 

y 

(pc) 


Rt 
(pc) 


nil — niu 
(Mq) 


(3 


m 
(Me) 


mi5 
(Mq) 


(Myr) 


101 


2 


4 


2 


X 


10'* 


30 


1.11 


1-150 


0.526 


5.04 


84.3 


2.7 


102 


2 


4 


2 


X 


10^ 


30 


1.11 


0.1-150 


0.722 


0.73 


72.1 


2.9 


103 


2 


4 


2 


X 


10* 


30 


1.11 


1-75 


0.659 


4.37 


54.9 


4.0 


111 


2 


1 


2 


X 


10^ 


30 


1.11 


1-150 


0.526 


5.04 


84.3 


2.7 


112 


2 


7 


2 


X 


10^ 


30 


1.11 


1-150 


0.526 


5.04 


84.3 


2.7 


113 


1.5 


4 


2 


X 


10* 


30 


1.11 


1-150 


0.266 


12.25 


111.6 


2.3 


114 


2.5 


4 


2 


X 


10^ 


30 


1.11 


1-150 


0.763 


2.76 


49.5 


5.1 


115 


1.5 


4 


2 


X 


10^ 


30 


1.11 


0.1-150 


0.430 


3.87 


109.9 


2.2 


116 


2.5 


4 


2 


X 


10* 


30 


1.11 


0.1-150 


1.000 


0.29 


26.7 


11.2 


117 


2 


1 


2 


X 


10^ 


30 


1.11 


0.1-150 


0.722 


0.73 


72.1 


2.8 


118 


2 


7 


2 


X 


10^ 


30 


1.11 


0.1-150 


0.722 


0.73 


72.1 


2.8 


121 


1.5 


4 


1 


X 


10^ 


30 


1.89 


1-150 


0.966 


12.25 


140.3 


3.0 


122 


2 


4 


1 


X 


10^ 


30 


1.89 


1-150 


1.020 


5.04 


127.3 


4.1 


123 


2.5 


4 


1 


X 


10^ 


30 


1.89 


1-150 


1.124 


2.76 


92.9 


15.8 


124"* 


1.5 


4 


5 


X 


10^ 


30 


0.70 


1-150 




12.25 


62.4 


1.2 


125 


2 


4 


5 


X 


10^ 


30 


0.70 


1-150 


0.108 


5.04 


42.7 


1.6 


126 


2.5 


4 


5 


X 


10^ 


30 


0.70 


1-150 


0.384 


2.76 


24.2 


2.4 


127 


2 


4 


1 


X 


10^ 


30 


1.89 


0.1-150 


1.032 


0.73 


119.5 


4.8 


128 


2.5 
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0.29 
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0.73 
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0.29 
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5.04 
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1.1 
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134 


2 


4 


5 


X 


10^ 


100 


1.43 


1-150 


0.108 


5.04 


42.7 


4.1 
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4.9 
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3.89 
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5 
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0.70 
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0.306 
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2.2 
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3.19 


59.4 


4.0 
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1 


X 
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30 


1.89 
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1.083 


3.19 


105.6 


8.6 
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5 


X 


10^ 


30 


0.70 


0.1-150 


0.735 


0.36 


16.9 


3.3 
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4 


2 


X 


10^ 


30 


1.11 


0.1-150 


0.913 


0.36 


37.6 


6.4 
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0.36 
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40.8 
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3.89 


1-150 


1.083 


3.19 


105.6 


24.2 



becomes a negative value for this model, so an equal logarithmic binning was used instead. 
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Table 2 

Properties of Compact Young Clusters 



Cluster 


Log M 


(pc) 


Log pav 

(M0PC-3) 


Age 
(Myr) 


Arches 


4.3 


0.2 


5.8 


1-2 


Quintuplet 


3.8 


1.0 


3.2 


3-5 


R136 


4.5 


1.6 


3.3 


< 1-2 



Note. From Tabic 5 of Figcr ct al. (1999). M is the 
total cluster mass in all stars extrapolated from obser- 
vation down to a lower-mass cutoff of 1 Mq, assuming a 
Salpctcr IMF slope and an upper mass cutoff of 120 M©. 
Rav is the average projected separation from the controid 
position. Since the low-mass end has not been identified 
in these clusters yet, the radius may be larger than the 
value given here when mass segregation is present, pav 
is M divided by the volume inside Rav 
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Fig. 1. — The effective lifetime Tg// of a star with ini- 
tial mass to(0) defined by equation (^. Schaller et al. 
(1992) was adopted for m{t). A logarithmic extrapo- 
lation has been made for m(0) > 120 M©. These Teff 
values are used for Tj^±i in equation (^. 



Fig. 2. — Evolution of M and pc for models 111 (Wq = 
1; dotted hues), 101 (Wq = 4; dashed lines), and 112 
(Wo — 7; solid lines), which have the same initial 
conditions except for Wq. These models have lagely 
different tec but very similar tev This shows that t^v 
is independent of Wq. 
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-3 -2 -1 



Log R / R^ 

Fig. 3. — Density profiles of mass component 1 to 15 
(bottom to top) of model 142 at f = 2 Myr (postcol- 
lapse phase). Both core and envelope are dominated 
by heavy components. 
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Fig. 4. — The ratios of tgv between models with the 
same parameters except mi {tev for m; = 0.1 Mq over 
tev for mi = IMq). Models shown here have Rg = 
30 pc and Wq = 4. Triangles are for M = 5 x 10^ Mq, 
squares for M = 2 x 10'* Mq, and pentagons for M = 
10^ Mq. The dotted line represents t^v ratios of unity. 
For a < 2, tev is nearly independent of mi in range 
0.1 < mi/M0 < 1. For M < 2 X 10* Mq, tev differs 
only by a factor of < 2 up to a = 2.5 in the same m; 
range. 
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Fig. 5. — tev as a function of mi5 (a) and a (b) for models with mi — 1 Mq and Rg ~ 30 pc, and Wo=4. Symbols 
are the same as in Figure |[ Dashed line represents tev oc nii^ relation. 
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Fig. 6. — tev as a function of Rg for models with mi — IM0, Wq — 4, and a = 2 (a) and 2.35 (b). Symbols are 
the same as in Figure ^. Without stellar evolution, tev oc (dashed lines) would be expected. 
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Fig. 7. — pc (solid lines) and M (dashed lines) evo- 
lution of model 142 (thin lines) and the same model 
with a ten times larger Rt (thick lines). The density 
profile of the latter is set to be the same as the former, 
so the cluster of the latter model does not initially fill 
up the tidal radius. This figure clearly shows the ef- 
fects of strong tidal fields on t^v ■ The pc of the model 
with a ten times larger Rt has a second core collapse 
at the end of the evolution. 



> 0.5 




o 



O 



Fig. 8. — Evolution of pc (solid line), M (dashed 
line), Rt (dotted line), and Vc (dash-dotted line) of 
model 142. M and Rt are normalized with their ini- 
tial values. The estimated locations in t-pav plane 
for the Arches and Quintuplet clusters are indicated 
with bars (data from Table |^; note that p^s of the 
clusters will be slightly higher than pav expecially for 
the Arches which shows a high central concentration). 
The Arches, the more compact of the two, is situated 
near the epoch of core collapse, and the Quintuplet, 
which looks much less bound, corresponds to the fi- 
nal disruption phase of model 142. The Figure may 
imply that two clusters had similar initial conditions, 
and represent two distinctly different epochs in their 
evolution. 
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Fig. 9. Evolution of a for the whole cluster (thick 
hne) and four inner-to-outer, equally-spaced annuli 
between the center and Rt (thin lines; from bottom 
to top) of model 142. Considering the lower detection 
limit in observations and stellar evolution of massive 
stars, a was obtained for the mass range of [3, 30] M©. 
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